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Abstract 

We study a two component dark matter candidate inspired by the Minimal Walking 
Technicolor model. Dark matter consists of a dominant SIMP-like dark atom component 
made of bound states between primordial helium nuclei and a doubly charged technilep- 
ton, and a small WIMP-like component made of another dark atom bound state between 
a doubly charged technibaryon and a technilepton. This scenario is consistent with di¬ 
rect search experimental findings because the dominant SIMP component interacts too 
strongly to reach the depths of current detectors with sufficient energy to recoil and the 
WIMP-like component is too small to cause significant amount of events. In this context 
a metastable technibaryon that decays to and can in principle explain 

the observed positron excess by AMS-02 and PAMELA, while being consistent with the 
photon flux observed by FERMI/LAT. We scan the parameters of the model and we find 
the best possible fit to the latest experimental data. We find that there is a small range 
of parameter space that this scenario can be realised under certain conditions regarding 
the cosmic ray propagation and the final state radiation. This range of parameters fall 
inside the region where the current run of LHC can probe, and therefore it will soon be 
possible to either verify or exclude conclusively this model of dark matter. 

Preprint: CP3-Origins-2015-033 DNRF90, DIAS-2015-33. 


1 Introduction 

Dark matter in the form of dark atoms can be a viable solution to the current experimental 
puzzle. Even stable electrically charged particles can exist hidden in such atoms, bound by 
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ordinary Coulomb interactions (see [DEIE] and references therein). Stable particles with 
charge -1 are excluded due to overproduction of anomalous isotopes of hydrogen. However, 
for negatively doubly charged particles the constraints are much weaker and this possibility 
remains open. 

There are several models where heavy stable -2 charged species O can emerge naturally: 

(a) AC-leptons, from theories beyond the Standard Model, based on the approach of almost- 
commutative geometry in El El E] • 

(b) Technileptons and anti-technibaryons in the framework of Walking Technicolor (WTC) 

in m nni EH En Ea EH. 

(c) stable “heavy quark clusters” IJUU formed by anti-t/ quark of a 4th generation IIESI 

EH EH EH EH]. 

All these models predict also corresponding +2 charge particles. If these positively charged 
particles remain free in the early Universe, they can recombine with ordinary electrons in 
anomalous helium, which is strongly constrained in terrestrial matter. Therefore a cosmological 
scenario should provide a mechanism which suppresses anomalous helium. There are three 
possible mechanisms that can provide a suppression: 

(i) The abundance of anomalous helium in the Galaxy may be signihcant, but in terrestrial 
matter a recombination mechanism could suppress it below experimental upper limits 
HE). This possibility requires a new U(l) gauge symmetry which materialises as a long 
range interaction among dark matter particles mediated by a dark photon. 

(ii) Free positively charged particles are already suppressed in the early Universe and the 
abundance of anomalous helium in the Galaxy is negligible PEE]. 

(hi) Excessive positively charged particles are created in dehcit relative to excessive negatively 
charged particles and bind with them in neutral atom-like states 1201EU. 

These three possibilities correspond to three different cosmological scenarios of dark atoms. The 
hrst one is realized in the scenario with AG leptons, forming neutral AG atoms [5]. The second 
assumes a charge asymmetry of the O which forms the atom-like states with primordial 
helium P EE]. In the third case, which can be realized in WTG, the excess of positively 
charged technibarions UU~^^ is much smaller, than the excess of negatively charged ( , so 

that all the free UU~^~^ are captured by ( in (^ UU^~^) neutral bound states [201 El]- 

In WTG (see [S] and references therein for a review) new stable species (techniquarks and 
technileptons) that transform under nontrivial representations of the SU(2) electroweak group, 
participate in sphaleron processes that eventually determine the relative abundance between 
baryons and -2 charged species [H EH EH EH EH 1211 EH [2E] . 

After formation during the Standard Big Bang Nucleosynthesis (BBN) era, screens the 
O charged particles in composite ) OHe “atoms” [TEj. In all the models of OHe, 

O behaves either as a lepton or as a specihc “heavy quark cluster” with strongly suppressed 
hadronic interactions. However OHe is a strongly interactive massive particle (SIMP) due to 
the ^He nucleus. The cosmological scenario of the OHe Universe involves only one parameter 
of new physics i.e. the mass of O . Such a scenario is insensitive to the properties of O 
(except its mass), since the main features of the OHe dark atoms are determined by their nuclear 
interacting helium shell. In terrestrial matter such dark matter species are slowed down and 
cannot cause signihcant nuclear recoil in the underground detectors, making them elusive in 
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direct WIMP search experiments (where detection is based on nuclear recoil) such as CDMS, 
XENONIOO and LUX. The positive results of DAMA and CoGeNT experiments (see [27j for 
review and references) can End in this scenario a nontrivial explanation due to a low energy 
radiative capture of OHe by intermediate mass nuclei [IllZlE]. This can explain the negative 
results of the XENONIOO and LUX experiments. The rate of this capture is proportional to the 
temperature: this leads to a suppression of this effect in cryogenic detectors, such as CDMS. 

OHe collisions in the central part of the Galaxy can potentially lead to OHe excitations 
and de-excitations with pair production in EO transitions that can explain the excess of the 
positron-annihilation line, observed by INTEGRAL in the galactic bulge [21 [3l [2H |28] . 

One should note that the nuclear physics of OHe is in the course of development. The basic 
feature that determines the elastic OHe-nuclei interactions is an electric dipole potential due to 
the OHe’s charged constituents. This potential requires a specihc quantum mechanical analy¬ 
sis [29] . The lack of a Coulomb barrier seems to lead inevitably to overproduction of anomalous 
isotopes [30]. However one should study precisely how quantum tunneling is modified in this 
case. 

It has been shown 0 mi IMl [23 ES] that a two-component dark atom scenario is also 
possible and can be naturally realized in the framework of WTC, where both stable doubly 
charged technileptons ( , playing the role of O , and positively doubly charged technibaryons 

UU are present. The two-component dark matter scenario is realised by having an abundance 
of ( , which binds with and a much smaller excess of that binds with ( , 

forming a WIMP-like UU( component of dark matter. In the framework of WTC, UU~^~^ can 
be metastable, with a dominant decay channel to a pair of positively charged leptons. We 
argued in izu that even a 10 ® fraction of such positively charged techniparticles with a mass 
of ~ 1 TeV and a lifetime of ~ 10^° s, decaying to and can potentially 

provide an explanation for the observed excess of cosmic ray positrons, while being compatible 
with the observed gamma-ray background. 

The anomalous excess of high-energy positrons in cosmic rays was hrst observed by PAMELA 
EU and was later conhrmed by AMS-02 [32]. These results generated widespread interest, since 
the corresponding effect cannot be explained by positrons of only secondary origin and requires 
primary positron sources, e.g. annihilations or decays of dark matter particles. Recently 
the AMS-02 collaboration has reported new results on positron and electron fluxes in cosmic 
rays [33] and positron fraction [3l]. These measurements cover the energy range 0.5 — 700 GeV 
for electrons and 0.5 — 500 GeV for positrons and provides important information on the origins 
of primary positrons in the cosmic rays. In particular, the new results show for the first time, 
that the positron fraction no longer exhibits an increase with energy above ~ 200 GeV. The 
possibility of explaining the cosmic positron excess by decays of UU particles, comprising the 
tiny WIMP-component of dark matter in the considered scenario, was discussed in detail in [2T] . 
Here we make a blind search within the whole parameter space and we hnd the optimal values 
that give the best possible £t to both AMS-02 data on cosmic positron flux and FERMI/LAT 
data on diffuse gamma-ray flux [35] - 

The paper is organised as follows: in section we briefly discuss the high-energy positron 
abundance in cosmic rays and the calculation of cosmic positron and gamma-ray fluxes from 
decays of dark matter particles and in section we search for the best-fit parameter values 
and discuss the prospects of positron anomaly explanation in the framework of the considering 
model. We conclude in section HI 
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2 Positron and Photon Flux due to UU decays 

The abundance of high-energy positrons in cosmic rays (also known as “positron anomaly”), 
which was hrst observed by HEAT |3S] and was later conhrmed by various experiments, such as 
ATIC [57j, AMS-01 |3B], Fermi [32], PAMELA [31] and AMS-02 [32], can be a possible astro- 
physical manifestation of metastable dark matter. Alternative explanations of this phenomenon 
include old supernova remnants m ans pulsars m- Here we are going to consider the case 
where an unstable, though long-living, doubly-charged heavy particle which is part of a 

tiny ( 10 WIMP-like component of dark matter UU(, decays as UU —)■ ^ 

in principle with different relative probabilities. All decay modes produce positrons and pho¬ 
tons either directly or through decay cascades. Photons appear in leptonic decays as hnal state 
radiation (FSR), arising from r decays to hadrons and from inverse Compton scattering (ICS) 
on the interstellar electromagnetic background. The latter processes are not taken into account 
in this work, since their contribution was found to be subdominant [21] • Since there is no 
production of antiprotons in our scenario, diffuse gamma background gives the most stringent 
constraints. 


Here we shall briefly remind the calculations of positron and gamma-ray flux, done in [21] . 
The positron flux can be estimated a^ 

E)),E,. ( 1 ) 

Here c is the speed of light, r is the lifetime of UU, riioc = ^ ■ (0.3 GeV/cm^)m^^ is the local 
number density of UU particles with ^ = 10“® being the fraction of dark matter in the form of 
UU(, muu is the mass of UU, m is the mass of the positron, dN/dEo is the number of positrons 
produced in a single decay (obtained using Pythia 6.4 [13]), /3 ~ 10“^® s“^GeV“^ is the average 
energy loss rate and the factor 




2A3 


3A3 


( 2 ) 


takes into account effectively the finite size of the magnetic halo {h ~ 4kpc, R ~ ISkpc), A 
is the characteristic diffusion length, rj denotes a Heaviside step-function (see [21] for details). 
Here we will be interested in energies > 20 GeV, where the effects of solar modulation are 
negligible [33]. The spectrum of secondary positron was taken from [33]. 

For gamma rays produced by UU decays in our Galaxy, the flux arriving at the Earth is 
given by 

= —(3) 
T dvrAHobs iAOob, ^loc dE 

where we use an isothermal number density distribution = (skpc) +(8^kpc) ^ / g^^d r 

are respectively the distances from the Earth and the Galactic center. We obtain the averaged 
flux over the solid angle AGobs corresponding to Galactic lattitudes |6| > 20° and Galactic 
longitude 0 < / < 360°. Out of our Galaxy decays of UU are homogeneously distributed over 

^An accurate comparison to the semi-analytical approach, presented in [42) . didn’t point out any significant 
difference. 
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the Universe and should also contribute to the observed gamma-ray flux. This contribution is 
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where x = 1100 corresponds to the recombination epoch, (nmod) is the current cosmological 
number density of Uf/, ~ ln~^ inverse value of the Hubble 

parameter with t^oA being the age of the Universe. Ha and = 1 — Ha are respectively the 
current vacuum and matter relative densities. The total gamma-ray flux arriving at the Earth 

is 

E.,{E) = E^^\E) + E\^\E). (5) 

In [21] we showed that the considered model can in principle accommodate nicely the AMS-02 
data on positron abundance in cosmic rays [32], being consistent with the FERMI/LAT data 
on diffuse gamma-ray background (inclusive photon spectrum) [33] • In the following section we 
And the best fit parameters, taking into account the latest relevant data from AMS-02 [SS] |3l] 
and FERMI/LAT [15]. 


3 Fitting AMS-02 and FERMI Data 

In order to And the best possible fit for the data we vary four parameters: the mass of the 
decaying particle M, it’s lifetime r and two values of branching ratios Br^e and (the 

third one is automatically fixed by Bvrr = 1 — Br^e — Br^^). To And the parameter values 
which provide the best fit of the positron flux data and do not contradict the diffuse gamma- 
background data, we minimize the following function 


X = 


Nd.o.f. 


'N„, 


E 

2 = 1 


ams 


Fpos{E^^g, a) 
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Nf. 




i=i 
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- ^ermi) 
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( 6 ) 


Here {Eams,^ams) is a set of AMS-02 (energy, positron flux) data points {Nams = 37), 
Epos{E,d) denotes the predicted positron flux (see Eq. [^, aams is the AMS-02 statistical error, 
{Efermi^ ^fermi) is a Set of FERMI/LAT IGRB data points {N fermi = 7 - we choose only data 
points in the high-energy region, where the FERMI/LAT data give the most stringent limit 
on the predicted gamma-ray flux), and E^{E,d) stands for the predicted gamma-ray flux (see 
Eq. [^. O'fermi correspouds to the FERMI/LAT statistical error, and 9 (E^{E fermi, o) — ^ fermi) 
is a Heaviside step function (it excludes the contribution of the points where E^ < ^fermi)- 
To simplify this analysis we evaluate the minimal value of for each value of mass M in 
the range 500 — 1200 GeV (with a 100 GeV step). Once the mass is fixed, the parameters we 
scan over are a = {t, Bvee, Br^^p}. The number of statistical degrees of freedom is Nd.o.f. = 
Nams + Nfermi — 3 — 1. The lower limit on the mass value comes from the search for long-lived 
doubly charged particles at the LHG [16]. The rest of the model parameters can have any 
reasonable value (branching ratios should have a positive value up to 1 and the lifetime should 
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Figure 1: as a function of particle mass. For each value of the 1/1/ mass we have hxed the 

residuals parameters to the values that provide the best £t of the data (see text). Horizon¬ 
tal lines show the p-value for the corresponding x^, which can be roughly considered as the 
probability of the viability of the scenario. 



Figure 2: The positron fraction in the cosmic rays from decays of dark matter particles 
(red curve), corresponding to the best-£t parameters (M = 700 GeV, r = 8 ■ 10^°s, Hrge = 
0.182, = 0.394, Bvrr = 0.424), and the fraction of secondary positrons (gray line), com¬ 

pared to the latest AMS-02 data [3l] (blue dots). 
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Diffuse gamma background 



Figure 3: Gamma-ray flux multiplied by from decays of dark matter particles 
in the Galaxy and beyond (green curve), corresponding to the best-£t parameters 
(M = 700 GeV, r = 8 ■ 10^° s, Br^e = 0.182, = 0.394, Bvrr = 0.424), compared to the lat¬ 

est FERMI/LAT data on isotropic diffuse gamma-ray background |15] (|&| > 20°, 0° < / < 360° 
with point sources removed and without diffuse emission attributed to the interactions of Galac¬ 
tic cosmic rays with gas and radiation helds (foreground)); here three different foreground mod¬ 
els A (red dots), B (blue dots) and G (yellow dots) are shown). In our analysis we have used 
model B. 

be around 10^° s). By using the aforementioned method we obtain for each value of UU mass 
the set of optimal parameters that £t the observational data in the best way. To probe whether 
or not a particular scenario can provide a good agreement with the experimental data, we treat 

as a common statistical test and use the Pearson’s chi-squared test. 

The results of our analysis are presented in Fig. As it can be seen the best possible 
scenario takes place at = 1.236 for M = 700 GeV, r = 8 ■ 10^°s, Rree = 0.182, = 

0.394, Bttt = 0.424. However, it is clear that the obtained parameters provide a not-so-good 
description of experimental data due to the constraints, following from FERMI/LAT data. For 
this set of parameters, we show respectively in Fig. and the predicted positron fraction and 
gamma-ray background, compared to the experimental data. 

Despite the fact that the fit is not good, there are ways that can alter the picture and 
render the considered scenario viable. One possibility is that FSR production can turn out to 
be suppressed in the case of two identical particle final state (like e+e"*"), compared to the case, 
where two non-identical particles (like e'''e“) are created, as it was assumed in our calculations. 
A second way out is related to the current uncertainties of the cosmic-ray propagation model, 
that includes errors of diffusion coefficients, unaccounted effects such as re-acceleration and 
choice of source distribution and halo geometry. A rough estimate shows (see Fig. that 
each of the two proposed ideas applied to the considered model, can significantly reduce the 
y^ value for a considerable range of UU masses. However, a more detailed study and better 
understanding of the uncertainties of cosmic rays propagation is required in order to clarify the 
viability of the considered UU decay model. In addition, one should note that inverse Gompton 
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as a function of mass 



Figure 4: as a function of UU mass for the two possibilities described in the text: the green 

(lower) curve is obtained under the assumption that the theoretical values of the flux have a 
20% error in all the points (due to the uncertainties in cosmic-ray propagation); the magenta 
(upper) curve corresponds to the case where FSR is turned off (i.e. the case of two identical 
particles in the hnal state, where FSR could be signihcantly suppressed). 


scattering processes were not taken into account and their contribution into diffuse gamma-ray 
flux can also affect slightly the £t. 


4 Conclusions 

If the Standard Model can serve as a guiding theory for the explanation of the dark matter 
puzzle, there is no reason, why dark matter must be a single particle [2]. Multi-component dark 
matter with components that contribute different amounts to the total dark matter abundance 
is a possible scenario. There can easily be a dark matter realization where a dominant SIMP- 
like component like He^ is accompanied by a small WIMP-like component UU(. This can 
be particularly interesting because there have been several issues with the validity of the Cold 
Collisionless Dark Matter paradigm recently. The core vs cusp problem in dwarf galaxies, the 
number of satellite galaxies and the “too big to fail problem” might be manifestations of the 
fact that dark matter might be warmer than what was thought until now (see e.g. 03 and ref¬ 
erences therein). A two component dark matter scenario can in fact resolve the aforementioned 
issues. The scenario considered here involves a form of neutral OHe dark atoms made of stable 
heavy doubly charged particles and primordial He nuclei bound by ordinary Coulomb interac¬ 
tions. This scenario can be realized in the framework of Minimal Walking Technicolor, where 
an exact relation between the dark matter density and baryon asymmetry can be naturally 
obtained. Strict conservation of technilepton charge together with approximate conservation 
of technibaryon charge results in the prediction of two types of doubly charged species with 
strongly unequal excess, i.e. dominant negatively charged technileptons ( and a strongly 

















sub dominant component of technibaryons , bound with ^ in a sparse component of 

WIMP-like dark atoms {( UU^~^). Direct searches for WIMPs put severe constraints on the 

presence of this component. However we have demonstrated in |2T] that the existence of a 
metastable positively doubly charged techniparticle, forming this tiny subdominant WIMP-like 
dark atom component and satisfying the direct WIMP searches constraints, can play an im¬ 
portant role in the indirect effects of dark matter. We found that decays of such positively 
charged constituents of WIMP-like dark atoms to leptons e'^e'^, can explain the 

observed excess of high energy cosmic ray positrons, while being compatible with the observed 
gamma-ray background. These decays are naturally facilitated by GUT scale interactions. The 
best £t of the data takes place for a mass of this doubly charged particle around 1 TeV or 
below making it accessible in the current run of LHC. Our rehned best £t analysis of the recent 
AMS-02 and FERMI/LAT data presented here, can provide a crucial test for the decaying dark 
atom hypothesis especially since LHC in its second run will be able to probe the parameter 
space where this dark atom model is viable. 
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